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Abstr act. We study integrability properties of a general version of the Boltzmann 
collision operator for hard and soft potentials in n-dimensions. A reformulation of 
the collisional integrals allows us to write the weak form of the collision operator as 
a weighted convolution, where the weight is given by an operator invariant under 
rotations. Using a symmetrization technique in L p we prove a Young's inequality 
for hard potentials, which is sharp for Maxwell molecules in the L 2 case. Fur- 
ther, we find a new Hardy-Littlewood-Sobolev type of inequality for Boltzmann 
collision integrals with soft potentials. The same method extends to radially sym- 
metric, non-increasing potentials that lie in some L^ eak or L s . The method we 
use resembles a Brascamp, Lieb and Luttinger approach for multilinear weighted 
convolution inequalities and follows a weak formulation setting. Consequently, it is 
closely connected to the classical analysis of Young and Hardy-Littlewood-Sobolev 
inequalities. In all cases, the inequality constants are explicitly given by formulas 
depending on integrability conditions of the angular cross section (in the spirit 
of Grad cut-off). As an additional application of the technique we also obtain 
estimates with exponential weights for hard potentials in both conservative and 
dissipativc interactions. 

1. Introduction 

1.1. Background. The nonlinear Boltzmann equation is a classical model for a gas 
at low or moderate densities. The gas in a spatial domain ft C R™, n > 2, is modeled 
by the evolution of the mass density function f(x, v, t), (x, v) <E fl x R n , modeling the 
probability of finding a particle at position x, with velocity v at the time t € 1. The 
transport equation for / reads 

(d t + v-V x )f = Q(fJ), (1.1) 

where Q(f, f) is a quadratic integral operator, expressing the change of / due to instan- 
taneous binary collisions of particles. The precise form of Q(f, f) will be introduced 
below, for both conservative (elastic) [20] and dissipative (inelastic) interactions [19] . 
The Q(f, f) operator factorizes as the difference of two positive operators, usually 
denoted by the Q + (f,f)(x,v,t) rate of gain of probability due to two pre-collisional 
velocities for which one of them will take the direction v and the Q~(f, f)(x, v, t) rate 
of loss of probability due to particles that get knocked out of the direction v. 

In addition, these operators depend on the form of their collision kernels which 
model the collision frequency depending on the intramolecular potentials between in- 
teracting particles. More specifically, these kernels depend on functions of the relative 
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speed and on the scattering angle, the latter modeled by an angular function referred 
as the angular cross section. In all the cases we assume that the angular cross section 
is modeled by an integrable angular function on the S'™ -1 sphere (this condition, in 
the theory of the Boltzmann equation, is called the Grad cut-off assumption). The 
collisional kernels are further divided into the following classes: hard potentials, corre- 
sponding to unbounded forms of the relative speed, modeling stronger collision rates, 
and soft potentials modeling weaker collision rates, both as the relative speed is larger; 
and Maxwell molecule type of interactions where collisional kernels are independent of 
the relative speed. 

1.2. Aim of the paper. It is the purpose of this work to investigate the L r -integrability 
(in velocity) of the gain operator as a bilinear form Q + ( f, g)(v) acting on probability 
mass densities f(v) and g(v), and to search for exact representation formulas for the 
inequality constants and possible optimal estimates depending on the L p and L q norms 
of / and <?, respectively. This provides a new detailed study of the Boltzmann equa- 
tion with harmonic and functional analysis tools. It is possible due to the weighted 
convolution nature of its collisional integral operator in weak form. The work focus on 
questions on sharp constants for estimates in such functional spaces as well as on the 
existence and exact description of maximizers. 

In order to achieve our goals, we introduce a representation of collisional integrals in 
weak form that allows us to write the gain term of the collision operator as a weighted 
convolution, with the weight being given by a suitable bilinear operator invariant under 
rotations, acting on the test function. 

Following the initial idea developed in [3] , the main ingredient to approach the con- 
volution estimates is an L p -radial symmetrization technique, a genuinely new addition 
to the Boltzmann theory, used here with much more generality. The authors showed 
the Young's inequality for the elastic hard potential case and calculated exact con- 
stants which were proved optimal for Maxwellian molecule type models in L 2 . The 
method was applied to the strong formulation of the gain operator in the Carleman 
representation. 

The new representation for the weak formulation of the gain collisional integral as 
a bilinear form Q + {f,g) as a symmetric weighted convolution presented here allow us, 
very handily, to extend their results to collisional forms with soft (i.e. singular) and 
symmetric decaying potentials, as well as to the case of dissipative interactions. We 
also observe that the corresponding loss bilinear form Q~(f,g) does not have the same 
convolution symmetry property as the gain one, and so the results we obtain reflects 
this discrepancy between the gain and loss part of the Boltzmann collision operator. 

More specifically, we employ still here the L p -radial symmetrization technique, but 
in addition we develop a Brascamp, Lieb and Luttinger-type approach to multilinear 
weighted convolution inequalities associated to a weak formulation of the collisional 
forms [TTl [18] . In addition, like in the connections between the works of Beckner [7] , 
Brascamp-Lieb [El [28] and multilinear convex inequalities such as the sharp Young's 
and the sharp Hardy-Littlewood-Sobolev's, we also obtain Young's inequality for hard 
potentials, and the Hardy-Littlewood-Sobolev for soft potentials where the weighted 
convolution structure has a singular kernel. We also show a Young's inequality for ra- 
dially symmetric, non-increasing potentials, where the weighted convolution structure 
contains a singular kernel. In our case, the invariant- under- rot at ions weight multi- 
plier in the convolutional structure of the weak formulation of the collisional operator 
introduces a non-linear change of coordinates due to the angular integration in the 
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(n — l)-dimensional sphere. That means our inequalities and methods are not a di- 
rect application of those in [TTl [TBI [28] but rather an analog result for dissipative 
Boltzmann operators of collisional type. The range of convexity exponents depends 
on the convolution symmetry property of the associated bilinear form, and recovers 
the Hardy-Littlewood-Sobolev convexity relation in the gain operator Q + case, while 
is restricted for the loss operator Q~. 

All these estimates are valid for conservative or dissipative interactions between par- 
ticles, whose restitution coefficients have absolute continuity and monotonicity proper- 
ties as functions of the impact parameter (see (j2.4[) below). 

Finally we obtain Young's inequalities with exponential weights for hard potentials, 
both for the elastic and strictly (dissipative) inelastic case. It is remarkable that, in the 
dissipative e(z) < 1 case, tails of order A and decay rate a are preserved (called stretched 
exponential tails), with order < A < 2. However, we show that in the elastic case 
these estimates hold for classical Maxwellian tails but they have a polynomial weighted 
norm in one of the components. These last estimates are of interest for the study, for 
example, of L , L°° and tail propagation dissipative kinetic collisional models, such as 
granular flows. 

In all cases, the inequalities constants are given by explicit formulas depending only 
on certain integrability conditions of the angular cross section (in the spirit of Grad 
cut-off) . 

In short, we hope that this work, inspired from the harmonic analysis toolbox, 
provides simplified methods, extensions to the full range of exponents including soft 
potentials and radially non-increasing potentials, and also extensions to dissipative 
interactions and weighted estimates, to achieve a better qualitative understanding of 
all the convolution-like inequalities governing the Boltzmann theory. 

2. Main Results 

2.1. Preliminaries. In this paper we study the integrability properties of the Boltz- 
mann collision operator in the case of elastic or inelastic collisions. The gain part of 
this bilinear collision operator, commonly denoted by Q + (f,g), is non-local in both 
components. It is defined via duality by 

/ Q + (f,g)(vMv)dv:= [ [ f(v)g(v*) f ^(v')B(\u\,u ■ u) du dv, dv, (2.1) 

where the functions f,g 7 ^p £ Co(K") (continuous with compact support). The symbol 
u represents the unitary vector in the direction of u (u — u/\u\) and duj is the surface 
measure on the sphere S" 1-1 . The variables (pre-collision velocities), v',v'^ (post- 
collision velocities) and u (relative velocity) are related by 

u = v — v* , v = v — — (tt — \u\lj) and v + v t =v' + v'^ (2.2) 

The corresponding loss part of Boltzmann collision operator, is only non-local in one 
of its its components, say g. It is defined in a strong form by 

Q-(/,$)(«)HI&IUs»-i)/(«) / g(v*Mu)dv*. (2.3) 

We notice that, wile Q + is a bilinear operator with non-locality in both its components, 
Q~ is non-local only in its second component g and maintains locality in its first 
component /. This difference in the nature of there non- locality plays a crucial role 
in the range of exponents for the Hardy-Littlewood-Sobolev inequalities for the Q~ 
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with singular potentials (Corollaries l9l and fT0|) . The inelastic properties of the collision 
operator are encoded in the positive scalar function (3 : [0,oo) — > [5, 1] defined by 
f3{z) :— , where parameter e is the so-called restitution coefficient which enjoys 

the following two properties that assure micro-reversibility of the interactions: 

(i) z i — v e{z) is absolutely continuous and non-increasing. 

(ii) z i — v ze(z) is non-decreasing. 

The dependence of the restitution coefficient on the physical variables is commonly 



given by z = \u\ y ^ > i- e - the restitution coefficient e, and thus /3, depends only on 
the impact velocity 



13 f\ u \ = ^— <- . (2.4) 

We point out that the model for e could be more complex (for example assuming de- 
pendence on macroscopic variables like temperature), however this will not be the case 
in this paper. 

The particle interaction is elastic when the parameter (3 = 1 , and is referred as sticky 
particles when f3 = 1/2. A complete discussion of the physical aspects of the restitu- 
tion coefficient can be found in [19] . Standard models for the restitution coefficient, 
for example constant restitution coefficient and viscoelastic hard spheres, satisfy the 
assumptions (i) and (ii) above. We refer the interested reader to [2J, [5], [TB], [25] and 
[30] for additional numerical and mathematical references that use this class of models. 

The nature of the interactions modeled by Q + is encoded in the kernel B(\u\,u ■ uj) 
modeled by strength of intramolecular potentials, and many physical models accept 
the representation (henceforth assumed) 

-B(ju|, u ■ uj) — \u\ x b(u ■ uj) with — n < A. 

Depending on the parameter A the interaction receives different names: soft-potentials 
when — n < A < 0, meaning that larger relative velocity corresponds to a weaker 
collision frequency; Maxwell molecules type of interactions when A = 0, of collision 
frequency independent of the relative velocity; hard potentials when A > 0, mean- 
ing that larger relative velocity corresponds to stronger collision frequency. For the 
(nonnegative) angular kernel b(u ■ uj) we will require the Grad cut-off assumption 
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b{u ■ ui)Au) < oo. 

We refer to [16] and [23] for a detailed discussion on the inelastic collision operator. 



2.2. Description of the results. Very recently, Alonso and Carneiro [3] revisited the 
L p -analysis of the operator Q + in the elastic case (restitution coefficient e = 1) for the 
case of Maxwell type of interactions and hard potentials (i.e. < A < 1) and developed, 
by means of radial symmetrization, a Young's inequality approach that produces sharp 
constants. Their approach involved the analysis of the collisional integral on the gain 
operator in strong form by means of the Carleman integral representation. 

For our current goals, we now use a weak formulation of the collisional integral writ- 
ten as weighted convolution, where the weight is a nonlinear invariant-under-rotations 
operator acting on the test function (see (|2.8j) below). 

From now on we work in the setting of dissipative interactions satisfying conditions 
(i), (ii) and (|2.4j) as described above. 
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Let ip and 4> be bounded and continuous functions. Define the bilinear operator 

V(ip,4>)(u):= I i/j(u-)<p(u + )b(u ■ u) du) , (2.5) 
•As™- 1 

where the symbols u + and u~ , commonly known as Bobylev's variables, are defined by 

u~ := f(w - \u\lj) and u + := u - u~ = (1 - 0)u + |(u + \u\oS). (2.6) 

We shall see in Section 3 that, under conditions (|2.4[) and (|2.6[) , the operator (|2.5[) 
has a certain invariance under the group of rotations in R". This operator was first 
introduced by Bobylev [9] [10] in a slightly different setting where it was shown that 

Q+iL9)=V(f,g), (2.7) 

in the elastic Maxwell molecules case (i.e. A = and (3 = 1). More recently, when 
introducing the study of the Boltzmann dissipative model for Maxwell type of interac- 
tions in [T2], the authors showed that such relation also holds for constant /3 ^ 1. For 
references on the use of the Fourier transform in the analysis of the elastic Boltzmann 
collision operator one can consult |21] and |27] . In the dissipative interaction case we 
also refer to [25] for the Fourier representation of the Q + operator. 

In particular, from (|2.ip and (|2.5p we obtain the following relation between the 
operators Q + and V appearing in the weak formulation of the gain operator, now 
written in velocity and relative velocity coordinates 

(f,g)(v)i ) {v)dv= I [ f{v)g(v-u)V{T v Ki>,l)(u)\u\ x dudv 

Jr™ Jr™ ^2 g\ 

f(u + v)g(v)V(l,T- v ip)(u) \u\ x dudv, 



where r and 1Z are the translation and reflection operators 

T v ip{x) := ip(x — v) and TZip(x) := ip(—x). 

Representation (|2.8[) exhibits the nature of the weighted convolution structure of the 
weak formulation of the gain operator Q + and also shows that the integrability proper- 
ties of the collision operator Q + are closely related to those of the bilinear operator V . 
A similar approach was carried out in |24j which relates the operator Q + to a slightly 
different angular averaging operator. 

One of the advantages of working with the bilinear operator V is that a simple change 
of variables in (|2.8I) allow us to change the roles of / and g without essentially changing 
the convolution structure of the integral. This convolution symmetry will be specially 
important later on in the proof of the Hardy-Littlewood-Sobolev (HLS) inequality 
for Q + in the full range of exponents (Theorem [2]) . Such convolution symmetry is 
not present in the loss operator Q~ as we shall see in Corollaries l9l and ITOl that the 
corresponding weak formulation for Q~(f,g) has a convolution structure only in its 
second component g, but not in its first one /. As a consequence, the HLS inequality 
for this operator will only be valid in a restricted range (Section 6.2). 

In Section 3 we develop the L p -analysis of the operator V by means of the weak rep- 
resentation (|2.8[) using a radial symmetrization method introduced in [3 for a different 
representation of the gain operator in strong form by means of the Carleman represen- 
tation. This approach using (12.8[) . provides exact constants in some of our inequalities, 
both in the hard and soft potentials settings, that are sharp in some cases. In general, 
they are calculated depending on (explicit) integral conditions on the angular cross 
section. 
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The results of Section 3 are crucial for all results proved in the following three 
sections. 

In Section 4 we revisit Young's inequality for hard potentials, where the new addi- 
tions are the calculation of exact constants for conservative or dissipative interactions, 
which are proved to be sharp in the constant energy dissipation rate for Maxwell type 
of interactions and selected L p -exponents. 

To this end, consider the weighted Lebesgue spaces L^!(R n ) (p > 1, k > 0) defined 
by the norm 

II/IUjcb-) = (f IfW (i + M pfe )d« 

We prove the following. 

Theorem 1. Let 1 < p, q, r < oo with l/p + 1/q = 1 + 1/r. Assume that 

B(\u\,u-cj) = \u\ x b(u-Lj), 

with A > 0. For a > 0, the bilinear operator Q + extends to a bounded operato^ from 
L p a+x (R n ) x L q a+x {W l ) -> L r a (R n ) via the estimate 

\\Q + (f,g)\\ KiRn) <C H/ll ^ tR ,.) \\ghl +x (m)- (2.9) 
The constant C is given by 

r 

C = 2 A+2+ 2+7 \S n - 2 \ \ 2V J (i=£) 2r' d£(s)j 

r' 

(/jm+p-Aiwr 17 ^))^ 

n~ 3 

where the measure on 1] is defined as d£„(s) = 6(s)(l— s ) 2 ds and /3q = /3(0). 
In the case (p,q,r) = (1, 1, 1) the constant C is to be understood as 

a r 1 
C = 2 A+3+ 2 \S n - 2 \ J d^(s). 

In Section 5 we prove a Hardy-Littlewood-Sobolev-type inequality for the collision 
operator in the case of soft potentials, as follows: 

Theorem 2. Let 1 < p, q, r < oo with —n < A < and l/p + 1/q = 1 + X/n + 1/r. 

For the kernel 

B(\u\,u ■ uj) — \u\ x b{u ■ u>), 

the bilinear operator Q + extends to a bounded operator from LP (W 1 ) x L q (W 1 ) — s> Z7(R") 
via the estimate 

\\Q + (f,g)\\ Lr(Rn) <c\\f\\ LP(m \\g\\ Lq{Rn) . (2.io) 

The constant C is explicit in (|5.10|) , (|5.14j) and (|5.18[) . 



i In this paper we shall prove all the inequalities for a dense subspace of smooth functions giving the 
desired extensions for L p , 1 < p < oo. The L°° bounds, when applicable, are easily treated directly 
just pulling out the L°°-norm out of the integrals. 
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The constants we obtain for the two inequalities above are explicit, but generally 
not sharp. Only in the cases a = A = 0, (p, q, r) — (2, 1, 2) and (p, q, r) — (1, 2, 2) we 
find the sharp constant for the Young's inequality (|2.9[) (see Corollary [5]). In fact, the 
quest for the sharp forms of these inequalities in the other cases, which could be seen 
as analogues of the remarkable works of Beckner [7] , Brascamp-Lieb [17] and Lieb [28] , 
seems to be a very difficult problem in harmonic analysis. 

Finally, in Section 6, we provide further applications of the methods described here. 
We first present in Section 6.1 a description of the estimates for collision kernels with 
radial non-increasing potentials that leads to multiple Young's inequalities in Corollar- 
ies[7]and[8] In Section 6.2, we study the corresponding inequalities for the loss operator 
Q~ in Corollaries l9l and fTOl and see that the lack of convolution symmetry restricts the 
range of exponents. 

In Section 6.3 we prove Young-type estimates with exponential weights for hard 
potentials both for the elastic and strictly (dissipative) inelastic case. We show in 
Proposition[TT]that in the elastic case the weighted estimates with classical Maxwellian 
tails hold but have a polynomial weighted norm in one of the components. However, 
things get better in the strict dissipative e(z) < 1 case, as shown in Theorem 112] where 
tails of order A and decay rate a are preserved (called stretched exponential tails) , with 
order < A < 2. These are important tools in the study of propagation of moments 
[IB] and L 1 — L°°-exponentially weighted comparison principles [2"3"] . 

2.3. Related literature. Young-type inequalities in Theorem [T] for the collisional in- 
tegrals, in the case of Maxwell type and hard potentials, reveal the convolution nature 
of the operator Q + . In the elastic case, this observation was first introduced by Gustafs- 
son [26] for an assumed angular cross section function with pointwise cut-off from away 
from zero (grazing collisions) and 7r (head-on collisions). Later, Mouhot and Villani 
[31] Theorem 2.1] revisited the work of Gustafsson under the same conditions to obtain 
a different control of the constants (observe however that the constant in their Young's 
inequality blows up at the endpoints, due to the pointwise cut-off assumption). Also, 
standard integrability of the angular cross section is required away from the cut-off. 
Later, Duduchava-Kirsch-Rjasanow [22] used elementary techniques to obtain a simpler 
proof of the control of the L p -norm for the Q + operator in three dimensions by means 
of L 1 and L p norms where the constants in our present work match precisely the ones 
obtained in [52] for the case (p, q, r) = (p,l,p). About the same time, Gamba, Panferov 
and Villani, (24] Lemma 4.1], presented related estimates that use dual integrability 
of the functions, namely, / and g belonging to L 1 n L p . This stronger assumption 
presented in [24] allowed them to remove the restriction of pointwise cut-offs and use 
only the integrability of the cross section. Related work on L\ estimates for the elastic 
or the dissipative collisional integrals was also done by Gamba-Panferov- Villani [2"3"] . 
Bobylev-Gamba-Panferov |16) . and Mischler-Mouhot-Ricard |30) . 

In particular, the following observation is worth stressing: some techniques valid in 
the elastic case are not available in the inelastic (dissipative) one. For instance, defining 
the angular kernel in a half spherical domain is not possible due to the lack of symmetry, 
and in particular, such lack of symmetry forces a dual integrability assumption in 
the derivation of L p estimates done in [24, Lemma 4.1]. In fact, dual integrability 
assumption is adequate for the quadratic operator, and also in some cases for bilinear 
estimates, such as when studying propagation of derivatives for the conservative case. 
However, it is not always the case that such stronger assumptions hold, an so one needs 
to study a purely bilinear estimate. This is precisely the case for comparison techniques, 
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where a function g, which in principle is unrelated to a solution /, is "compared" with 
the solution and so the bilinear form Q{g, f) appears. In general, there is no reason 
on why this function g may have dual integrability, and so, the previous technique will 
not apply. Finally, we mention that the authors in |16j have already used the angular 
averaging mechanism occurring for V (to be introduced in the next section) in order to 
obtain sharp moment decay formulas for the gain operator, both in the case of elastic 
and inelastic (with constant restitution coefficient) hard potentials. These results lead 
to the study of the regularity and the aforementioned Gaussian propagation for the 
solutions of the Boltzmann equation using comparison techniques [1] and |23j . 

In the soft potentials case, Theorem [2] also reinforces the convolution character of 
Q + {f, g)- It is however important to notice that the weak formulation of the collisional 
integral as a weighted convolution (representation (|2.8p ) is crucial to attain the results 
by estimating multilinear integrals by convex type estimates, maximized in their sym- 
metrization and the connections to Brascamp-Lieb-Luttinger inequalities [T5| |T7] . It 
is also worth mentioning that the classical Hardy-Littlewood-Sobolev was used in con- 
nection to Cancellation Lemmas of Alexandre-Desvillettes-Villani-Wennberg [1] and 
Villani |32j who constructed estimates with singular kernels and non-integrable cross 
sections, assuming L p k integrability to control the i^-norm of the total collisional oper- 
ator. Our work is not related to theirs, but rather extends to any radial non-increasing 
potential, showing that the collisional integral satisfies a Hardy-Littlewood-Sobolev- 
type inequality both for Q + and Q~ independently (see also Corollaries l9l and [T0|) . 
providing exact representations for the constants, and assuming only L p -integrability. 
In this sense it is the first time where such connection was made. 

Finally, recent main applications of this work are mentioned below. Alonso and 
Gamba [5] used this result to obtain classical solutions and instability , 1 < p < oo, 
for the Cauchy problem associated to the Boltzmann equation for soft potentials and 
singular radially symmetric decreasing potentials, with integrable cross section, for 
initial data near vacuum or large data near local Maxwellian distribution. In addition, 
Alonso and Lods [6] have shown the pointwise control to inelastic homogeneous cooling 
solutions for the inelastic Boltzmann equation for hard spheres and the conjectured 
Haff Law by means of a corresponding Young-type estimate for the inelastic Boltzmann 
collision operator with stretched exponential weights. 



Let G = SO(n) be the group of rotations of R" (orthonormal transformations of 
determinant 1), in which we will use the variable R to designate a generic rotation. We 
assume that the Haar measure d/i of this compact topological group is normalized so 
that 



3. Radial Symmetrization and the Operator V 




Let / <E L p (W l ), p > 1. We define the radial symmetrization /* by 



l 




(3.1) 



and 



fto(x) =ess sup lyl=lxl \ f(y)\ 



(3.2) 
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where the essential sup in (|3.2p is taken over the sphere of radius |x| with respect to the 
surface measure over this sphere. This radial rearrangement /* defined in (|3.1j) - (|3.2[) 
can be seen as an L p -average of / over all the rotations R G G and it satisfies the 
following properties: 

(i) /* is radial. 

(ii) If / is continuous (or compactly supported) then /* is also continuous (or 
compactly supported). 

(iii) If g is a radial function then (fg)*(x) — f*(x)g(x). 

(iv) Let dv be a rotationally invariant measure on R™. Then 



\f(x)\*>dv(x) = / \fi(x)\*d V (x) 



In particular, 



LP( 



WfZ\ 



We noticed that this is not the classical rearrangement by monotonicity of level sets 
associated to /, stills has the property that preserves the L p -norm, and so we are 
driven to prove a convex inequality type of result that can be viewed as an analog to 
a Brascamp Lieb and Luttinger type inequalities, where multilinear convolution type 
integrals are maximized on their (classical) rearrangements and controlled by multiple 
Minkowski type of convexity estimates [18j [17] . 

The following result, our first in this manuscript, is a redo of [3] Lemma 4]. We 
produce here a different proof of that in [3] (which use a Carleman integral type repre- 
sentation) , and present a new approach that exhibits the similarity to those arguments 
of classical rearrangement inequalities such as those in [18j El] • These new approach is 
crucial to perform the extensions to soft potentials through Hardy-Littlewood-Sobolev 
and Lieb [29] type of inequalities. 

Lemma 3. Let f,g,ip€ Cq(W) and 1/p + 1/q + 1/r = 1, with 1 < p,q,r < oo. Then 



V(f,g)(u)ip(u)du 



< 



r(f;,g*)(u)r r (u)du. 



Proof. From (12.4[) . (|2.5[) and (|2.6I) we observe that for any rotation R one has 

V(f 1 g){Ru)=V(foR,goR){u). 

Therefore, 



/ 



V(f,g)(u)i>(u)du 



V(f,g)(Ru)^(Ru)du 
T{foR,goR)(u)ip{Ru) du 



(3.3) 



< 



gn-l 



\f{Ru~)\ \g(Ru + )\ \tp{Ru)\b(u-uj)dwdu. 



Note that the left hand side of (|3 . 3[) is independent of R. Thus, an integration over the 
group G = SO(n) leads to 



V(f,g)(uMu)du 



< 



\f(Ru-)\ \g(Ru + )\ \ip(Ru)\dfi(R) b(u-uj)dLudu 



G 



(3.4) 
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An application of Holder's inequality with exponents p, q and r yields 

\f(Ru-)\ \g(Ru+)\ \1>(Ru)\dii(R) < #(u~) g*(u + ) <(«), 

which together with equation (|3.4I) proves the lemma. □ 

Lemma [3] shows that L p -estimates for the operator V will follow by considering 
radial functions. If / : R™ — > R is radial, we define the function / : R + — > R by 

/(*) = /(M). 

In addition, for any p > 1 and a £ R we have 

/•OO 

/(ir) p \x\ a dx = jS*™" 1 ! / f(t) p t n - 1+a dt. (3.5) 
Jo 

Hence, if we define the measure v a on R™ by 

di/ Q (x) = |x| Q dx , 

and the measure er" on R + by 

do%(t) =t n - 1+a dt, 

equation (|3.5I) translates to 

l 

||/||iP(E",dj/ a ) = ||/||l,P(K+,dO- (3-6) 

In the following computation we show how the operator V simplifies to a 1-dimensional 
operator when applied to radial functions. If / and g are radial, then 

V(f,g)(u) = [ f(\u~\)g(\u + \) b(u-u)du; 
jS" -1 

/(ai(|u|,u- uj))g(a 2 (\u\,u- u))) b{u ■ u)duj (3.7) 

s „-i 

/l „_3 
^ /(ai(|«|, a ))5(o2(|«|, s))b(s) (1 - ,s 2 )— ds. 

The functions a\ and a 2 are defined on R + x [—1,1] — > R + by 

ai (x,s)=(3x(^) 1/2 and a 2 (z, s) = x [(1±2) + (1 - /3) 2 (±^)] 1/2 . (3.8) 
We conclude from (13.71) that 



Flf.gH.r) = |S"' -| I f{a x {x,s))~g{a 2 {x,s)) d#(«) > (3.9) 



where the measure ^ on [—1,1] is defined as 



n— 3 



d^(s) = 6( S )(l-^)— ds. 

In virtue of equation (|3.9p we define the following bilinear operator for any two bounded 
and continuous functions /, g : R+ — > R, 

B(f,g)(x):= J f( ai (x,s)) g(a 2 (x,s))d£(s). (3.10) 

Remark. It is worth to notice that in the case of constant parameter (3 (which includes 
clastic interactions) the functions a± and a 2 of the variable interactions are actually 
functions of the form a\ — xai(s) and a 2 = xa 2 (s); that is a\ and a 2 are first order 
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homogeneity in their radial part and their angular part is a positive, bounded by unity 
function of the angular parametrization s. This observation is at the heart of the 
analysis made in |14j which shows compactness properties of the spectral structure 
associated to the bilinear form (|3.10[) . 

For the operator in f|3 . 10[) we have the following bound. 

Lemma 4. Let 1 < p,q,r < oo with l/p + l/q = 1/r. For f E L p (R + ,da°) and 
g e L q (R+,do-%) we have 

P(/)P)IL-(R+,d CT «) ^ C \\f\\LP(M+,da S ) lblli»(a+,dcr«) > ( 3A1 ) 

where the constant C is given in (|3.16|) below. In the case of constant restitution 
coefficient e, corresponding to a constant parameter f3 = (1 + e)/2, one can show that 

n+q /■! n+a n+a 

C(n,a,p,q,b,p)=p- P J (^)~ 2 *> [(^) + (1 - Pf ~ ^ d&(«) 

(3.12) 

is sharp. 

Proof. Using Minkowski's inequality and Holder's inequality with exponents p/r and 
q/r we obtain 

1 

\f( ai (x,s))\ r \g(a 2 (x,s))\ r daZ(x)Y dtiis) 



W,g)\ 



L''(R+,dcr°) 
1 



< 



< 



-1 



\f{ ai (x,8))\ P doZ(x) 



15(03(0;, do^^) de„(.s) 



Since the function z — > ze(z) is non-decreasing, the change of variables y — ai(x, s) is 
valid for any fixed s 6 [—1,1), and its inverse Jacobian satisfies 

da! 



dx 



Moreover, using the fact that /3 > 1/2, we arrive at 



n + Q 



n-\-a 



(3.13) 



(3.14) 



Using a similar analysis for the change of variables y = a2{x,s), exploiting the fact 
that j3 is non-increasing, we obtain 



dao 



dx 



> 



[(!^) + (l-/3 ) 2 (^)] 



(3.15) 



where /3q = f3(0). We then arrive at 



J^lgia^s^da^y < [(^) + (1 - A,) 2 (*?)] 



2g 



llfl'llL«(R+,do-«) ' 



This gives (|3.11j) with constant 

n+q /*! n+a n+a 

C = 2 v J 2 ? [{^) + (l-M 2 (^)V 2q (3-16) 

In the case of constant (3, the Jacobians (|3.13p and (|3.15p can be explicitly computed 
and the proposed change of variables leads to the constant (13. 12)) . To prove that the 
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constant (13.12[) is the best possible in this case, one can consider the sequences {f e } 
and {g e } with e > defined by 

e i/p x -( n + a -e)/p for o < x < 1 , 
otherwise. 



and 

9e(x) = 

Clearly, 



e l/q x -(n+ a -e)/q foj . Q < X < 1 

otherwise. 



H/e|lLP(K+,dO — ll5cllL9(K+,d<r°) — 1 ' 

and one can check that 

\\B(fe,9e)\\ L r(J l + iAa a ) -+C, 

as e — )• 0, where C is the constant defined in (|3.12[) . The detailed argument is outlined 
in [3], in the case /3 = 1. □ 

From Lemma [3] we have 

\\P{f,g)\W{M^,dv a ) < WP{fp,9q)\\L-(R".du a ) , 

where l/p + l/q = 1/r. Using equations (|3.6|) . (|3.9p and Lemma|l]we obtain 
\\V(r pl g* q )\\ L r {Rn , d „ a) = IS"- 1 !- Kg^) 



= |5n-x| P \\B(f*,g*)\\ L r (M+!d<) (317) 

i 

< C |S m ~ 1 | r jS" 1 " 2 ] ||/^||iP(M+,dO [|5glU9(R+,dff») 
= C \S n 2 | ||/||z,J>(Mn !<il/a ) ||si|z,<r(]R'»,dl' ct ) ) 

and thus we have proved the following result. 

Theorem 5. Let 1 < p, q,r < oo wii/i l/p+l/q = 1/r, and a£l. TTie bilinear oper- 
ator V extends to a bounded operator from L p (M. n ,dv a ) x L q (M. n ,dv a ) to L r (M. n , dv a ) 
via the estimate 

Ik (/>flOlll,»-(B n ,dWa) - C ll/II.L*"(B n ,di/ ,) IIS , Hz,«(E",dt'a) - 

Moreover, in the case of constant restitution coefficient e, the constant 
zs sharp. 

An application of Theorem [5] with Fourier transform methods provides sharp esti- 
mates for the (2,1,2) and (1, 2, 2) Young's inequalities in the case of Maxwell molecules 
and constant parameter (3 (they will be treated with more generality in the next sec- 
tion). These are the only cases where we are able to explicitly find the sharp constant 
and exhibit a family of maximizers for the Young's inequality. 

Corollary 6. Let f 6 L 1 (M n ) and g € L 2 (R n ). Then 

\\Q + (f,9)\\ L2{Rn) < Co U/IUip.) ||ff||i»(R»), (3.18) 
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with the sharp constant given by 



Co = IS" 



[(i±.)+(l_ 0)3 (!=£)]-* d£ (s) . 



Similarly, for f G L 2 (M n ) and g G L 1 (K n ) we have 

||Q + (/,ff)|| i 2 (R „ ) < C\ ||/||z,2(R«) ||5|Ul(R«), 

with the sharp constant given by 

a = \s n - 2 1/3-2 / (i^) 4 d ^ (s) . 



(3.19) 



V(f,9) 



L 2 (R") 



(3.20) 



Proof. To prove Q3.18P observe that 

\\Q + (f,g)\\ LHRn) = \\Q+(Lg) 

< Co ||/||l°°(R") llffllL 2 (R") < C ||/||ii(R») ||5||l 2 (R")) 

with the constant Co coming directly from Theorem [5] To guarantee that Co is indeed 
the sharp constant in the inequality (I3.20[) we need approximating sequences f e and 
g e slightly different from those presented in the end of the proof of Lemma [H since we 
would like to impose the additional constraint / > to have H/Hl^r™) = II/IIl^r™)- 
Heuristically, this can be done by considering / = S (x) the Dirac delta and so / = 1. 
In practice we should choose / £ a Gaussian approximation of the identity by putting 



AW 



and 



e l/2 r (n-e)/2 for q < x < 1 , 

otherwise. 



A similar consideration applies to the inequality (|3.19l) . 

4. Young's Type Inequality for Hard Potentials 



□ 



The goal of this section is to prove Theorem [T] First we treat the case a = A = 0. 
The main idea is to use the relation (|2.8I) that establishes a connection between the 
operators Q + and V, and then use the L p -knowledge of the operator V from the 
previous section. In what follows we assume that all the functions involved are non- 
negative (since when working with L p -norms one can always use the modulus of a 
function). From (12.81) we have 



/ := 



Q + (f,g)(vW(v)dv 



f(v)g(v -u)T(T v Kip,l){u)dudv. (4.1) 



Suppose first that (p,q,r) ^ (1, 1, 1), (1, oo, oo), (oo, 1, oo). The exponents p,q,r in 
Theorem[T]satisfy 1 jp' -f-l/g' + l/r = 1, and thus we can regroup the terms conveniently 
and use Holder's inequality 



I = 



f(v) P rg(v-U)r} (/(«) 1' V{T V U^, 1)(«) * 



/ jl zL 

(g(v-u)P'V(T v Kip,l)(u)p' 



(4.2) 



dudv < Iihh, 
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where 



1 

h := ( / / f(vfg(v -ufdudv' 



h~[f / f(v) p V(T v Tlip,l)(u) r dudv 

g{v - u) q V{T v Tl^, \)(u) r ' dudv 



(4.3) 



g(v) q V(l, T^){u) r dudv 

Recall that r and 1Z are unitary operators in the LP spaces, thus, from (|4.2[) and 
Theorem [5] we obtain 

I < C\\f\\ LP{Mn) \\g\\ L «(W.n) WWf'iRn), 

with constant given by 



f"l _r±_ \ q' 

r=\S"-'\[2- I (1=*) ar'd£*00 



(4.4) 



which concludes the proof. The cases we left over are easier. Indeed, for (p, q, r) — 
(1,1,1), it is a matter of pulling the L°°-norm of V(T v TZtp,l)(u) out of the integral 
(|4.1[) using Theorem [S] with (00,00,00), thus arriving at the constant 



C=\S n - 2 \ J 



In the case (p,q,r) — (00, 1, 00), it is just a matter of pulling out the L°°-norm of / 
out of the integral (|4.1[) . and repeating the process as in integral J 3 above. The case 
(p, q, r) — (1, 00, 00) is basically the same. In both these cases, the final constant can be 
read from (|4.4p . Note that one would obtain the same constants starting with the easy 
cases (p, q, r) = (1, 1, 1), (1, 00, 00), (00, 1, 00) and using Riesz-Thorin interpolation, but 
we chose to do this directly. 

In the case where a + A > and A > 0, we use two additional inequalities in order to 
control the post collisional local energy and powers of the post collisional velocity and 
of the the relative velocity by powers of the integration variables appearing on each 
= 1, 2, 3 in (|4.3[) . Indeed, from the energy dissipation we have 

KI 2 + KI 2 <M 2 + kl 2 , 

and thus 

| v - u-\ a = \v'\ a < (M 2 + H 2 )" 72 < 2 Q / 2 (\v\ a + \v- u\ a ) . (4.5) 

Also, 

\u\ x < (\v ~u\ + M) A < 2 A (\v - u\ x + \v\ x ) . (4.6) 
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Now, take as a test function ip a {v) = ip(v)\v\ a . The first step is to use (|2.8|) and (|4.6|) 
to obtain 



Q + (f,g)(v)ij a (v)dv= / f(v)g(v- u)V(T v H^ a , 1)(«) \u\ x dudv 

JR™ JR™ 

<2 A {/ / f(v)g(v-u)r(T v Kil} a ,l)(u) \v\ x dudv (4.7) 

- u)r{T v nip ai \){u) \v - u| A dudu| . 
Making use of (l2~5j) and (|4~5|) yields 



V(T v lZip a , l)(u) = / ip(v — u )\v — u \ a b{u-ui)duj 

<2 Q/2 |y ^(w-u^)|w| Q 6(u-w)do; (4.8) 

+ J ij](v — u~)\v — u\ a b(u • w)dw| . 
Recalling the notation ip a (v) = ^(v)!?^ , combining (|4.T[) and (|4.8p one obtains 

Q+(/, s )(«)i(.)d« 

<2 A+ f(/ y / A+0 («)ff(t;-u)^(T^,l)(«)d«dt) 

fx(v)g a (v - u)V(T v 1Zip, 1)(«) diidv (4.9) 
f a {v)g\(v - u)V(T v TZip, l)(u)dudv 



f(v)g\ +a (v - u)V{T v TZip, l)(u) dudvj 



We now repeat the procedure for the case a = A = 0, breaking each of the 4 integrals 
of (|4.9[) into 3 parts according to (|4.3p . One should then estimate each piece using 
Theorem [5] and simple inequalities of the form 

||/]|z,»(R«) < ||/|Us +A (K») and ||/ Q || L p(R") < ||/||lj +a (r«) ■ 
At the end we arrive at 

dv < 2 A+2 +! C ||/[[z£ +A (*») ll<?IU« +x(R ") Uhr'^y 
This proves that 

\\Q + (f,g)(v)\v\ a \\ Lr{&n) < 2 x + 2 +fc\\f\\ Ll+ ^ n) \\g\\ Ll+x(Rn) . 
A similar reasoning provides 

\\Q + (f,g)(v)\\ L r (mn) < 2 A+1 c||/|| iPQ+A(R „) lblU^ +x(R n) , 

and finally 

HQ+a^MlL^) < 2 A+2+ f +^ C ||/|U S+ , (R n) ||5lUs +A (R-) , (4-10) 
with C given in (|4.4[) . The proof of Theorem 1 is now completed. 
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5. Hardy-Littlewood-Sobolev Type Inequality for Soft Potentials 



In this section we study the collision operator for soft potentials and prove Theorem 
[U We divide this proof in three parts, making use of the convolution symmetry of 
Q = (/,s)- 



5.1. The case r < q. From (|2.8|) we have 



~(f,9)(v)ip(v) dv 



f(v)g(v - u)V(T v Kip, l)(u) \u\ x dudv 
r v TZg{u)P(T v niJj,l)(u) \u\ x du ] dv. (5.1) 



f(v) 



Applying Holder's inequality and then Theorem [5] to the inner integral of (|5.1|) . with 
(p, q, r) = (a, oo, a), we obtain 



T v Kg(u)V{T v K^,l)(u) \u\ x du< \\r(T v TZyj,l)\\ La{mndux) \\T v Tlg\\ 



< Ci ll^^|| L a (R „ idiyA) \\T V Kg\\ L a>p» idvx) 



l/o 



(M B *[u| A )(«) ' {\g\ a ' *\u\ x ){v) 



l/a' 



where l/a + l/a 1 = 1 (a to be chosen later) , and the constant C\ given by 



C x = \S r ' 



w+A 
2 a 



Therefore we arrive at 



1 n+A 



i<cx I f(v) (i^i°*h a )(«) 1/Q (i 5 r'*| W | A )(«) 



l/a' 



du. 



(5.2) 



(5.3) 



Applying Holder's inequality in (|5.3p with exponents 1/p + 1/6 + 1/c = 1 (6 and c to 
be chosen later) we arrive at 

l/a' 



(5.4) 



i<Cx Wfhn^) llH^Nlliv.^) ||lsl a '*M ; 

We now use the classical Hardy-Littlewood-Sobolev inequality to obtain 

llivr*H A Lva (K n)<^ iivii^ (R ") (s-s) 

and 



\g\ a * |«P 



I") 



(5.6) 



where 



a 1 X , a' 1 A 
1 + - = - and 1 + — = 

ban c e n 



The constants C2 and C3 (generally not sharp) are explicit in (29J p. 106]. Finally 
putting together (|5.5|) and (|5.6|) with (|5.4|) we arrive at 



y(R~) ll5lL«'e(Bn) ll^lli^(R»)' 



(5.7) 



To conclude the proof of the theorem it would suffice to have in (|5.7|) the relations 
a'e = q and ati = r'. Now it comes the moment to choose our variables. All the 



INEQUALITIES FOR THE BOLTZMANN OPERATOR 



17 



inequalities we used above will be well-posed if the following relations are satisfied 



1 1 



7 = 1, 1 < a < oo 



a a 
111 

- + - + - = 1, 1 < 6,c < oo 
p b c 

a I X 



1 



b d n' 

a' _ 1 X 

c e n 



b > a, 1 < d < oo 
c > a', 1 < e < oo 



(5.8) 



a e = q 
ad = r' 

The last two equations determine d and e in terms of a. The remaining linear system (in 
the variables 1/a, 1/a', 1/6 and 1/c) in undetermined because of the original relation 



1 + 1 

p q 



x 1 



One can check that the choice 



and 



111/ A 

h = ---( 1 + - 

o r a \ n 



1 1 1/A 



c q a 
with any 1 /a in the non-empty interval 



1 1 
max < — ,1 

I r' q(l + £ 



provides a solution for (I5.8[) . If we define 



1 

< — < mm ■ 



D 1 = d Cl /a Cl /a ' 



(5.9) 



(5.10) 



using expressions (|5.2I) , (|5.5I) , (|5.6I) and a choice of a given by (|5.9p , expression (|5.7I) is 
plainly equivalent to 

1 1 Q + (/,<?) 1 1 < A ll/IU^K") ||fflU«(R-), (5-11) 

finishing the proof in this case. 

5.2. The case r < p. Here we make use of the convolution symmetry of Q + given by 



/:= / Q + (f,g)(v)i;(v)dv= / f(u + v)g(v)T(l,T_ v iP)(u)\u\ x dudv 

9(v) ( I T- V f(u) V{l,r- v iP)(u) \u\ x du) dv. (5.12) 
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We now proceed exactly as in the proof of the first case, noticing that the roles of / 
and g (and thus of p and q) are interchanged. We will find 

/ < Ci C'l /a Cl /a ||flr||i,,(Rn) ||/|| L o' e(R „) HVIU-^CM")- 
where C4 (the analogous of C\ ) is given by Theorem [5] as 

/l ra+A 

while C5 and Cq (the analogues of C2 and C3) are given by classical HLS inequalities 
as in (|5.5p and (15.61) , for a choice of a now in the non-empty interval 



In the end, defining 



—-,1 — - — r— > < — < min < — - — - — r— , — 7 i • (5.13) 



D 2 = C i Cl /a Cl /a ' (5.14) 



we will have 

\\Q + (f,g)\\ L r m < d 2 \\f\\ L v m hWm^)- (5.15) 

5.3. The case r > max{p, q}. The remaining range will be covered by the multilinear 
Riesz-Thorin interpolation [33J Chapter 12, Theorem 3.3]. Recall that the triple (p, q, r) 
satisfies 1 < p, q, r < 00 and 

i + i=(l + ^U-- (5-16) 
p q \ n J r 

For fixed A and n, with — n < A < 0, there exist triples (pi,qi,r), with r < q%, and 
(j>2, Q2,r), with r < p2, such that: 

l + I=(l + ^ + i = l + l. (5.17) 
Pi 9i V nj r p 2 q 2 

Therefore, for the triple (pi,Qi, r) we have bound (|5 . 11 1) and for the triple (j>2, 92, r) we 
have bound (|5.15l) . From (|5 . 16[) and (|5 . 1 T[) there is a i, with < t < 1, such that 

/ 1 1 1 \ , 71 1 1\ /ill 
t + 1-t = -,-,- 

VPi 91 r / \P2 92 ry \p q r 

and we can apply the bilinear version of the Riesz-Thorin interpolation to get 

\\Q + (f,9)\\ Lrm < DlDi 1 -** \\f\\ LW hh^y (5.18) 
This concludes the proof. 



6. Some immediate applications 

The radial symmetrization technique and the key ideas in the proofs of the basic 
inequalities for Boltzmann (Young's and HLS) described here have a deeper range and 
may be applied to other important estimates in kinetic theory. In this section we 
discuss two additional examples. 
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6.1. Non-increasing Potentials. Assume in the following discussion that the colli- 
sion kernel has the form 

B(\u\,u-u) = $(u)b(u-u>), (6.1) 
for some radial non-increasing potential $ : W l R + . Let us retake the discussion of 
Section 3 by defining the measure z/$ on R n by 

df$(x) = $(x)dx , 

and the measure cr* on R + by 

do*(t) = <£(*) f^dt. 

We claim that under this convention, Lemma|4]and Theorem[5]are still valid with these 
measures replacing dv a and da" in the statement of these results. This follows from a 
simple observation in the proof of Lemma 0] when one performs the change of variables 

y = a 1 (x,s) = /3x(if^) 1/2 

1 /2 

Indeed, note that since $ is non-increasing and < (3 (^j 5 -) < 1, then 

$(x) < %). 

Hence, using that f3 > 1/2, the estimate (|3 . 14[) changes in this context to 

l 

(jf S ))rdtT*( a .)J <2*> 2 P\\f\\ LP(VL+ , d<) . 

Same observation is valid for the change of variables y = 012(2;, s), in the sense that 

1 

QT \g(a 2 (x,sWda%txf) * < [(^) + (1 - A,) 2 (^)]~^ IMI i9(R+ , d<) , 

where (3q — (3(0). We conclude that, for a radial non-increasing potential, Lemma U 
and Theorem [5] hold with the (non-sharp) constant 

C(n,p,?,6,/9) = 2-p J_{^)-^ [(^) + (l-M 2 (^)r~«d£( s ). 

An immediate consequence of this observation is that Theorem [5] holds for any 
collision kernel (|6.ip with a potential that is radial non-increasing and belongs to a 
weak Lebesgue space $ £ L^ eafc (R™). A careful reading on the proof of this theorem 
together with the discussion above leads to: 

Corollary 7 (weakly integrable potential for the gain). Let 1 < p,q,r,s < 00 with 
1/p + 1/q + 1/s = 1 + 1/r. Assume that $ £ L^ eafc (IR n ) is a radial non-increasing 
function. For the collision kernel (|6.1[) the bilinear operator Q + extends to a bounded 
operator from L p (M. n ) x L q (R n ) -)• L r (R n ) via the estimate 

||Q + (/:ff)|| i r (R „) < c II^IIl^^jr") ll/lliP(R") HsHi^R") - 

with C as in the proof of Theorem [H 

An interesting fact, that might be suitable for applications related to soft potentials, 
is that if we impose the stronger condition that our radial non-increasing potential 
belongs indeed to L S (R"), then we can include the endpoints in the range of exponents 
for our inequality. To see this, it is just a matter of following the proof of Theorem 
[5] and use the classical Young's inequality for convolution whenever the classical HLS 
inequality was used (passages (|5.5p and (|5.6[) ). One should obtain: 
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Corollary 8 (integrable potential for the gain). Let 1 < p,q,r,s < oo with l/p + 
1/q + l/s = 1 + 1/r. Assume that $ £ L s (W n ) is a radial non-increasing function. For 
the collision kernel (|6.1[) the bilinear operator Q + extends to a bounded operator from 
L p (R n ) x L q (R n ) -+ L r (WL n ) via the estimate 

I|Q + (/'.9)Hl-(R") ^ C W®\\l»(M") II/IIw(r») II5IL«(b»). 
with C as in the proof of Theorem 

6.2. Estimates for the loss operator Q~(f,g). This discussion leads to the appli- 
cation of the previous results to the loss operator Q~(f,g) defined in (|2.3[) in strong 
form. The following results assert that Theorem [2] and its Corollaries [7] and [8] are also 
valid for this operator, however in a restricted range of exponents due to the fact that 
this operator lacks the convolution symmetry of Q + given by (|2.8|) (see expression (|6.2|) 
in the proof below). We also show a counterexample that exhibits the lack of full range 
of exponents property that the gain operator Q + enjoys. 

Corollary 9 (weakly integrable potential for the loss). Let 1 < p,q,r,s < oo with 
1/p + 1/q + 1/s = 1 + 1/r and r < p. Assume that $ £ ^Lafc(^ B ) * s a ra dial non- 
increasing function. For the collision kernel (|6.ip the bilinear operator Q~ extends to 
a bounded operator from L p (R n ) x L q (W l ) —> L r (W l ) via the estimate 

\\Q~(f,g)\\ L r {Mn) < C\ la C\ la H&Hii^-i) H $ l!z4 eafc (R") ll/llip(R«) IIsIIl«(r«) • 
with C5, Cq and a as in the proof of Theorem^ (second case). 

Corollary 10 (integrable potential for the loss). Let 1 < p,q,r,s < 00 with l/p + 
1/q + l/s = 1 + 1/r and r < p. Assume that 4> £ L s (W l ) is a radial non-increasing 
function. For the collision kernel (|6.1[) the bilinear operator Q~ extends to a bounded 
operator from LP(« n ) x L*(R n ) -» L r (W l ) via the estimate 

\\Q~{f,g)\\ Lr{Rn) < II $ L S (k«) ll/IL^E") lbllL«(H») • 

Proof. These estimates can be proved noticing that for any test function ip one has 
Q-(f,g)(vMv)dv = \\b\\ LHs ^ 1) f g(v)( f f(u)^(u)^(v-u)du)dv 

JR™ VH» / 2) 

< l|6|lzi ( 5-i) / 9(v)W\ a * t>)(v)} 1/a [(\f\ a ' * *)(«)] 1/a 'd«, 

where 1/a+l/a' = 1. Now it is just a matter of following the proof of the second case of 
Theorem [21 Note that estimates (|5.5|) and (|5.6j) are achieved using the classical version 
of HLS inequality when $ is weakly integrable. Meanwhile, when $ is integrable they 
are achieved using the classical Young's inequality. It is precisely in this step that the 
range of exponents differ in the two corollaries. □ 

Remark: One should not expect Corollaries [S] and [TU] to hold outside the ranges 
above. For instance, in Corollary |H] we could examine the situation p = q = r = s = 2, 
taking the potential $(ir) = \x\~ n / 2 and 

{x) = { \A~ n/2 (In Mr 1 if \x\>e, 
\ otherwise. 

In this situation we would have for any v, 

/ g(v*)$(u)dv* = / g(x)$(v — x)dx 
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/ \x\- n/2 (\n\x\y 1 \x - v\' n/2 dx 

J{\x\>e} 
1 



> 



2 n/2 



\x\ "(log \x\) 1 dx = oc , 



{ \ x I >max{e, 1 1?| }} 

that makes Q~ (/, g)(v) to blow up pointwise in (|2.3|) . and thus the bound of Corollary 
IH1 would not hold. 

For Corollary 1101 we could even consider a simpler example taking $(a;) constant 
(s = oo), p = 1 and q = r = oo. From (|2.3I) it is easy to see that the bound in 
Corollary [TU1 cannot hold in this case. 

6.3. Inequalities with exponential weights. In this subsection we will present a 
Young-type estimate for the Boltzmann collision operator with stretched exponential 
weights. Our motivation for presenting such estimates is the study of tail propagation, 
for instance, in solutions of the elastic Boltzmann equation, in self similar solutions for 
the inelastic Boltzmann equation, or, in solutions of the inelastic Boltzmann adding 
heating sources. 

It is well known that the homogeneous elastic Boltzmann equation propagates so- 
lutions with Maxwellian tails. More specifically, it was proved by Bobylev [11] that 
such solutions f(t,v), in 3-dimensions and for hard spheres, satisfy the Z^-Maxwellian 
weighted estimate if initially so, that is, 



if 



fo(y) 



L 1 



< oo then 



sup 

t>o 



/(*, 



"I 



L 1 



< oo for some a < ao- 



Later, a program based in a comparison principle was introduced in [23], for the n- 
dimensional case with variable potentials and integrable differential cross section, to 
obtain exponentially weighted pointwise estimates 



sup 



/(t,«)e°W : 



< oo, 



under the assumption that /o satisfies similar L°°-Maxwellian integrability and not 
necessarily with same rate a. One of the key ingredients in this program is the inequality 
(see Lemma 5 in [23] ) 



Q + (e" 



/)(«)< 



■A 7 



< C 



with k = fc(A, n), with k = for the particular case of hard spheres in three dimension 
(i.e. A = 1 and n = 3). This estimate is a particular case of the Young's inequality with 
Maxwellian weights. The techniques used in [23], and the subsequent applications to 
propagation of Maxwellian tails for derivatives [J , were based on the strong formulation 
of the collisional integral by the Carleman representation. 

Further in the studies of inelastic Boltzmann equation, it was shown that solutions of 
the homogeneous inelastic Boltzmann for constant restitution coefficient behave com- 
pletely different from those of the elastic. Thus, the previous L°°-Maxwellian estimate 
does not hold for them. Indeed, as the time passes by, the gas cools down to complete 
rest. In other words, the density converges to a Dirac distribution at v = 0. This result 
was first rigorously proved for the Maxwell type of interaction models by Bobylev, Car- 
rillo and Gamba [12] . One may obtain such convergence by a dynamical rescaling that 
leaves the collisional integral invariant usually referred as a self-similar transformation. 
It was actually shown that for the case of any dissipative model of Maxwell type of 
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interactions, such transformation yield to a static frame where the distributions have 
power law for their high energy tails, p"3 j [T5 ] [14 ] . 

However this picture is very different for hard potentials, where stretched exponential 
tails must be expected. Indeed, Bobylev, Gamba and Panferov proved that any solution 
to the self-similar transformed static problem of the inelastic Boltzmann equation has 
solutions in L 1 with a weight that corresponds to a non Gaussian exponential high 
energy tail |T6] . Their studies also showed similar behavior to other inelastic collisional 
models under heating sources, where their stationary solutions are shown to have L 1 - 
exponential weighted estimates with weights depending on the restitution coefficient 
and the heating mechanism. Later, the time propagation of L 1 -cxponcntial tails by 
self-similarity, in the hard spheres case, was shown by Mischler, Mouhot and Ricart 

m- 

All these estimates, rigorously proved in the inelastic case for hard spheres and 
constant restitution coefficient, yield 



sup 

t>Q 



/ s (M)e aW 



< oo 



where f s is the self-similar profile. In addition, for hard potentials < A < 1 one 
expects using formal arguments, 



sup 

t>o 



f s (t,v)e aH 



< oo . 



(6.3) 



Consequently, a companion issue arises on the availability of exponentially weighted 
pointwise estimates for self-similar profiles f s associated to dissipative collisional equa- 
tions, that is 



sup 



fs(t,v)e a ^ 



< oo . 



(6.4) 



Passing from (|6.3j) to (|6.4j) will require an inequality of the type of (|6.5p . 



Another source for motivation in searching L p -exponentially weighted estimates can 
be found in a forthcoming work [B] where different issues of the self similar asymptotic 
of inelastic collisions for variable restitution coefficient are treated. 

Motivated by the above discussion, we present two Young-type estimates for the 
Boltzmann collision operator with stretched exponential weights which preserve the 
rate of the exponential for elastic and inelastic interactions. In order to shorten nota- 
tion, fix a > and define for 7 > the exponential weight and the Maxwellian weight 
respectively as 

M an {v) :=exp(-aH 7 ) and M a (v) := M a ,2(v). 

The first result of this section is an inequality which is ideal in the study of pointwise 
Maxwellian tails for hard potential in the elastic case. The proof that we give works 
for both elastic and inelastic cases. 

Proposition 11. Let 1 < p, q, r < 00 with l/p+ 1/q = 1 + l/r. Assume that 

B(\u\,u ■ lo) = \u\ x b(u ■ uj) , 

with A > 0. Then, for a > 0, 

\\Q + (f,g) Ma 1 \\ Lr t Rn) < C ||/ M~ 1 \\lp(m») lis M^WlKm™)- (6.5) 



INEQUALITIES FOR THE BOLTZMANN OPERATOR 23 

Proof. Using ([2~T]) and (|2?2]) we obtain 

I:= I Q+(f,g)(v){M-^)(v)dv 

f(v)g(v-u) I (M~ l tp) (v 1 ) \u\ x b{u-uj)dujdudv. (6.6) 

From the energy dissipation we have \v'\ 2 + \vl\ 2 < \v\ 2 + v*| 2 , and thus 

AC V) < M-\v) M-\v*) MM)- (6-7) 
Using (|6.7[) in (|6.6|) we obtain 

I< I I {M-\f){v){M- 1 g){v~u) f ^j(v') M a (v:)\u\ x b(u-uj)dududv. 
Jr" Jr" is™- 1 

(6.8) 



Recall that 



1— u-o? 



2 

Therefore, using the fact that (5 > 1/2, we have the simple pointwise estimate, 

h a <2 a (t^) a/2 ^^-'(^r 2 (K+--i A +ki A ). 

Further, v'^ = v* + u~ , then 

M a (v:){\v*+u-\ x + \v*\ x ) <C A , a (l + M A ). (6.9) 
Using ()6.9|) we have 

^{v')M a {v',)\u\ X b{u-L0)dL0 

(6.10) 

< C x , a (I + \v,\ x ) / i]j(v-u-)b{u-uj)dLj, 
where we have defined 

b (u ■ w) := (r=kjj &(«•")• 
Using (|6. 10|) in expression f|6 . 8|) we arrive at (recall that v* — v — u) 

We now have arrived at the same expression given in (|4.1[) . with /(u) changed by 
(.M" 1 /) (u) and g(u) changed by (M~ x g) (v) (l + |u| A ). Repeating the argument for 
the Young's inequality in Section 3 we will conclude that 

\\Q + (f,g) ■A^o 1 || i r( R n) < C C X:a ||/ M-^xsw \\g M^h^, (6.11) 

with C given by f|4.4[) with b in place of b and C A ,a from (I6.10p . This concludes the 
proof. □ 

If we restrict ourselves and impose "strict dissipation conditions" to the model we 
can do even better. Indeed, note that both norms in the right-hand side of estimate 
(|6.12p below are free from any extra weight coming from the potential. This feature 
is achieved by imposing additional conditions (meant as "strict dissipation") on the 
restitution coefficient e stated in the theorem. Such conditions are not stringent since 
they are satisfied by standard models such as viscoelastic spheres or constant restitution 
coefficient e < 1. 
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Theorem 12. Let 1 < p, q, r < oo with 1/p + 1/q = 1 + 1/r. Assume that 

B(\u\,u-u) = \u\ x b{u-uj), 
with < A < 2. Then, for non-increasing restitution coefficient such that e(z) < 1 for 

z e (o, oo), 

\Q + (L9) ^CaL^ ^ c Wf ^«,aIU'(«») lis *C,aIU«(r-)- ( 6 - 12 ) 

TTie constant C :— C(n, X,p, q,b, ft) is computed below in the proof. In the important 
case (p,q,r) — (oo, l,oo) this constant reduces to 



C = C(n, A) £ [(i±£) + (1 - /3(0)) 2 (if*)] n/2 6^) d«, 



(6.13) 



with 



6/9 («) 



i+hWf)iV /2 



&(*): 



1 - S 



= V^l - /?(^)) 2 + /? 2 (z) + 2(1 - /3(x))/3(x)s , for x = ^ 

Proof. Let us introduce the classical center of mass-relative velocity change of coordi- 
nates 

U = - "t, — and ?? = (1 - j9)tt + j8w. 



One can readily verify the following standard identities 



v' = U + ^ti and E:=\v\ 2 



\vJ 2 =2\U\ 2 



2 



Therefore, a direct calculation shows that 

Kl 2 = £(i±4^) , 

Thus, for any test function ip with unitary L r -norm, 



• , , 2 W \\u\ 
with £ := 1 " 1 < 1 . 



I := 



Q + (f,g)(v) [M-^) (v)dv 

f(v)g(v-u) f (m-^) (v')\u\ x b(u-cj) 



duj du dv. 



In order to estimate /, we split the integral into two regions of integration: A — {\u\ < 
1} and its complement. The integral in the region A is estimated by 



I A < 



(■Ma'xf) («) {M-^g) (v - u )J snl h ^ ■ w ) duj du dv - ( 6 - 14 ) 

Indeed, this estimate follows from local dissipation of energy estimates (notice here the 
condition < A < 2 is used), namely 

\v'\ x < (M 2 + KI 2 ) A/2 < E x ' 2 < \v\ x + M A = \v\ x + \v- u\\ 
Hence, proceeding as in the proof of Young's type inequality in section 4, it follows 



I A < Ci 



fM-\ 



Lp(R") 



9M-\ 



L<J(R") 



(6.15) 



where the constant C\ is explicitly computed by formula (|4.4j) . 
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The second integral is slightly more involved. In order to control the integral Ia<= 
we use the inelastic interactions law and conditions (|2.4p to our advantage. First, 
observe that i9 is a convex combination of two unitary vectors and since j3 > 1/2, the 
magnitude of *d increases as j3 gets closer to 1. Moreover, /3 is non-increasing, thus in 
A c = {\u\ > 1} one has 







< 



\M* ' w )l = vA 1 " + /? 2 (x) + 2(1 - Z?^))/^ • w, 



1— u-c 



with a; = 

Next, note that by assumptions (|2.4[) on the restitution coefficient e, the magnitude 
| dp | ^ 1 except for u • u = 1. Therefore, the magnitude of t/ is controlled by 



W? <E 



< J5 



2 



Also note that 

| U |A < 2 A/2 s A/2 ; 

and thus, for we can estimate the integral 1a" as follows, 

I M < 2 x ' 2 l f(v)g(v-u) [ £ A / 2 exp (aE x ' 2 (l±^ X ' 2 \ ^{v')b{u-Lu) 

Finally, noting that 

£ A / 2 exp ^ A / 2 (i±^) 

= exp (a E x ' 2 ) E x ' 2 exp ( -aE x ' 2 ( 1 



dui du dv. 



A/2 



A/2 



< a 1 sup {xe x \ ( 1 

a:>0 



1 + 10/3 I 

2 



A/2 



exp 



one can estimate Ia" by 

I A c < 

2 x ' 2 {a e)- 1 (M-lf) (v) (m"^) (v - u) j ^ _^{v') bp{u ■ us) dududv 

KC^iae^jfM-^ 
for any < A < 2, where we have defined 

A/2^ - 1 



Lvt 



9M-\ 



, (6.16) 



hi*) ■= i- 



2 



&(*). 



The constant C2 is defined again by (|4.4[) with fyg replacing 6. We obtain the final 
constant C by adding the constants obtained in (|6.15l) and (|6.16l) . namely, C = C\ +C2. 
Note that bp has a singularity at s = 1, which in most cases of interest is at least of 
first order. In the case (p,q,r) = (00, l,oo) the constant can be taken as showed in 
(I6T31) . □ 
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Remark: One wonders then, if estimate (|6.12|) is also valid in the elastic case, or 
perhaps, even the also useful (but weaker) estimate 
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